    * Table A.7: introduce an additional control variable: distance to ilorin
	matrix b53 = b50[1,1..6], 0.0, 0.0, 0.003, 0.26, 0.001, 0.0, 0.0, 0.0, 0.0,  0.001, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0
	
	set seed 101
	mixlogit d_  dist_dw_hf_2 dist_22   nu_8d  alt2-alt4  nu_8dsq distjtoil , group(choicesit) id(idid) cluster(q01_10) rand(doctor_tot  generator2009 beds p_i1 p_i2 p_i3 p_i4 ) from(b53, copy) 

